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Abstract 

In this work we present the results of a complete calculation of the 
77 — > 7r°7r° amplitude to leading order in the large iV approximation 
(N being the number of Goldstone bosons) up to order m\/F 2 . The 
amplitude turns to be proportional to that of -k + -k~ — > tt°it . In spite 
of the fact that this factorization property cannot hold in general (as 
it was recently pointed out by Morgan and Penington), it appears 
here since in the large N limit only the I = J = channel contributes 
to the 77 — > 7r 7r° reaction. Moreover it seems to be a reasonable 
approximation in this case since it is possible to reproduce, as a pre- 
diction, the experimental data starting from a one-parameter fit of the 
7T7T scattering data. 
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1 Introduction 



In the last years a great deal of work has been devoted to the 77 — > 7r°7r° 
reaction. The reasons for that are the experimental data available from 
the the Crystal Ball Collaboration [1] and that this process, being a pure 
finite one- loop effect, is a very good test for the modern Chiral Perturbation 
Theory (x-PT) [2] which successfully describes many low-energy hadronic 
reactions. However, for 77 7r°7r°, the one-loop prediction [3] disagrees 
with the experimental data even near threshold. This is due to the pion 
rescattering which is very strong in the I = J = channel. Nevertheless it 
has been shown that the use of dispersion relation techniques [4] can provide 
a good description of this reaction and reconcile data with xPT [5]. More 
recently it has been shown that a two loop-computation produces a cross- 
section that agrees rather well with the available data and with dispersion 
theoretical calculations even substantially above threshold [6]. 

Another completely different approach that has been proposed for this 
reaction by Im [7] is the so called large N expansion (N being the number 
of Nambu-Goldstone bosons). This method was first used for elastic pion 
scattering in [8] in the context of the Linear Sigma Model (LSM) and in [9] 
in the Non Linear Sigma Model (NLSM). Pion mass effects have also been 
included in [10] for the LSM and in [11] for the NLSM. In this last work it 
was shown how it is possible to fit the J = pion scattering phase-shifts up 
to very high energies. 

Nevertheless Im's work has been strongly critizised by Pennington and 
Morgan in [12] because some of his assumptions are only true for point inter- 
actions which is not the case of pions. Since Im's computation is not in fact 
a complete large N computation, these assumptions cannot be supported 
and thus we consider Pennington and Morgan criticisms appropriate. In this 
note we clarify this issue by making a complete large N computation of the 
77 — > 7r°7T reaction. Our main conclusion is that at the end the original Im's 
result is correct at leading order of the large N expansion and to first order 
in ml/F 2 so both Pennington and Morgan on one side as Im on the other, 
are simultaneously right. In addition we will show how the large iV approx- 
imation provides a very good description of the Crystal Ball experimental 
data. 
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2 The Chiral Lagrangian in Standard S 3 Co- 
ordinates 



As it is well known, the low-energy interactions of pions and photons can 
be described by a U(l) em gauged NLSM in the context of modern xPT . 
This NLSM is based in the coset space SU(2) L x SU(2) R / SU(2) L+R = 
0(4)/ 0(3) = S 3 = SU(2). In order to define the large N limit we will 
extend the coset space to 0(N + 1)/0(N) = S N and we will define the 
corresponding lagrangian as follows. First we start from the LSM 

£ = -d^dy + VNFm 2 a (1) 
2 

where 4> T = (ni, tt' 1: 7r 2 , 7r 2 , .., n^, n' k , ttn, a) which is 0(N + 1) invariant if it 
was not for the last piece that explicitly breaks this symmetry to O(N) (with 
N — 2k + 1). This contribution will provide the mass term for the pions and 
other interactions forced by the chiral symmetry. With 7r = kn the above 
lagrangian can be written as: 

C = \d,n a d^ a + \d^' a d^' a + ^tt^tt + l -d,ad»a + VNFm 2 a (2) 

with a = 1,2,.., A;, or, in terms of the complex charged fields (p a = (n a + 
in' a )/V2, as 

C = d^* a d^ a + \d^°d^° + \dpadPa + VNFm 2 a (3) 

The electromagnetic interactions are included by introducing the correspond- 
ing covariant derivative 

d^ a - D^ a = d^ a + ieA^ a (4) 

and adding the photon kinetic term 

2 2 

+ e 2 AV:/ + ieA^l^ + h.c. - \f^F^ (5) 

where F^ = d^A^ — d^A v . Now we obtain the NLSM describing the low- 
energy pion dynamics by setting T = NF 2 = f 2 to force the system to 
live in the S N sphere i.e. 

vr 2 + a 2 = NF 2 (6) 
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where n 2 = Ea=i^ a2 + K ,a2 ). Thus we have a = (NF 2 - tt 2 ) 1 / 2 
and the lagrangian of the U(l) em gauged NLSM becomes 

£ 9 = ^ d^ + -d,ndn +- _ _ ^ 




(7) 

with y? T = (v?i, ...,¥» fc ). 

This way to write the lagrangian describing the low-energy pion dynamics 
corresponds to a particular coordinate choice on the coset space S N but 
in principle any other parametrization connected with this one through an 
analytical change of coordinates on the coset manifold will yield the same 
S matrix elements (but not necessarily the same Green functions) [13]. In 
particular, for the N = 3 case a very popular coordinate choice consists in 
defining the following SU (2) = S 3 matrix: 

E = e—f^- (8) 

The change of coordinates <p l = 4>' 1 (uj), with i = 1,2,3 relating both parametriza- 
tions of the coset space is given by /^E = 4 + icjfr 1 with 4 = ^ f% — 2 

where <p 2 = X)f=i ^ with <p T = [0i, 2 , 03, 4>a\ — [^i, ^i, ^cb c]- The above 
lagrangian written in terms of these new coordinates reads 



f f 1 

C = ^-TrD^ (D^) + + i^Trm 2 (E + E+) - -F^F"" (9) 

where 

Z^E = d^ + ie[Q,E}A^ 

Q = diag (2/3, -1/3) (10) 

As it was commented above, both lagrangians yield the same S matrix ele- 
ments and the same observables like phase-shifts, cross sections, etc. How- 
ever, the Green functions could be different. Indeed the computations can 
be much easier using one set of coordinates than another since the Feynman 
rules and the diagrams contributing to some given process will be different 
in general. In particular, the computation of the one-loop 77 — > 7r°7r° ampli- 
tude is simpler using the standard S* 3 coordinates than the E (chiral) ones, 
which were used in [3] 
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3 The one-loop scattering amplitude 

The 77 — > 7T7T transition amplitude is defined as 

A ( 77 _> 7r o 7r O ) 7r+7r -) = ^ {h) e ^ {h) M % n o^- (11) 

where ei M , and fci, /C2 are the polarization vectors and the momenta of the 
two photons. 

The two one-loop diagrams for the process 77 — > 7r°7r° using standard S" 3 
coordinates are shown in Fig.l.a. Note that when this computation is done 
with chiral coordinates there are two more Feynman diagrams [3]. In our 
case, using dimensional regularization, the diagrams are given by: 



2e 2 



{NF 2 ) 



2 1 s "™. 



dt 



9iiv 



(NF 2 ) 2 
+ [(fci, A*) " 



[s - ml 
— (*2,i/)] 



dl- 



(i + hf-ml] [(£-k 2 ) 2 -ml 
(2£ + k 1 ) fi (2£-k 2 ) u 



(H 2 -m\) (£ + h)-ml 



k 2 f - ml 



(12) 



with 



74-e/ 



(27T) 



4-e 



(13) 



It is not difficult to see that the above Af® u and A/"^ correspond respec- 
tively to the terms M a ^ v + M c ^ v and + -M^„ of the second reference in 



[3]. The complete matrix element is given by 



2e 2 



(nf 2 y 



m. 



di- 



(t 2 -ml) (i + hY-ml (£-k 2 y-ml 



(14) 



which is proportional to the xPT lowest order n + n — > 7r°7r° amplitude given 
by (s - m 2 n )/f 2 . 

As it is well known, the integral above is finite and therefore no renor- 
malization is needed. This is because there is no four derivative term in the 
gauged NLSM contributing to 77 — > 7r°7r°. Therefore there is no any possible 
counterterm to absorb the divergences appearing in the one-loop computation 
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and thus the consistency of xPT renders this one-loop amplitude completely 
finite. The result is given by 



M 



2e 2 



(NF 2 



s — m, 



r ) (l67r 2 ) 



9nu 



k 2>1 k\ v 
hk 2 , 



mi 



1 + {\nQ n - my 

s 



1 + a/1 - 



(15) 



From this amplitude one can calculate the total cross-section which is found 
to be 



a (77 - 
where 



°) 



a 

8^' 



1- 



4ml 



s 



(T 7T 7T 



ttV) (16) 



(To (VT + TT- -> 7T 7T ) 

which is the lowest order 7r + 7r~ - 



2\ 2 



(s - m 2 ) 



327r(iVF 2 )' s 
7r°7r° scattering cross-section and 



(17) 



f(s) = 2 



In 2 



4m| 
s 









7T 2 




In 2 | 




s 





1 



7T 



(18) 



Thus, from the above equations one can realize that, as it was remarked 
in the second reference in [3], the lowest order 7r + 7r~ — > 7r°7r° amplitude and 
cross-section factorize the one-loop 77 — > 7r°7r° ones. 

From the diagrams in Fig.l. it is also possible to compute the reaction 
77 — > 7T + 7r~. However, in this case we have a tree level contribution (Fig.l.b) 
and, in addition, we have to consider at the one-loop level (Fig.l.c) the 
contribution coming from the pion wave-function renormalization (Fig.l.d). 
In the reference frame where the photon polarization vectors ei, e 2 satisfy 



h ■ e 1 = ki ■ e 2 = k 2 ■ e 1 = k 2 ■ e 2 
the amplitude at tree level is given by 







a tree 
(77^ 7r+7r ~) 



2ie < e 1 ■ e 2 ; — 

I k 1 -p+ 



V ■ • £2 
k 2 ■ p + 



(19) 



(20) 
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where fci,A; 2 are the photon momenta and p + ,p are the tt + , it outgoing 
momenta. Now including the one-loop corrections. 



a one— loop 

(77^7T+7T~ ) 



-ie 



lQn 2 NF 2 



{ei • e 2 (s + rn^ln 2 *^) 

p+ ■ e x p~ • e 2 p~ • tip + ■ e 2 



ei • e 2 - 



with 



A (m 2 ) = m 2 (A-lnm 2 ) 



A = - + log4vr-7 + l + ln/x 2 



so that the full amplitude is 



-loop 



^ (77 - vr+vr-) = [4^-) + ^^V.-). 
where the wave function renormalization constant is given by 

z - = 1 + mh^^ ( m 9 

so that we can write 

A ( + _\ . 2 / P + -eiP~-e2 P~-eiP + -e 2 

.A 1 77 — > 7r^7r J = lie < aei • e 2 — 



where a is given by 

a(s) = 1 - 



k\ ■ p A 



1 /s m 2 2 N 
16vr2iVF2 U + 2 V \ 



A; 2 • p A 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



4 The large N limit approximation 

In order to compute the 77 — > 7r°7r° cross-section at leading order in the 
1/iV expansion we go back to the lagrangian in eq.7. The pion scattering 
amplitudes to order 1/N and the lowest order in m 2 / / 2 were computed in 
[11]. The contributing diagrams in this approximation are shown in Fig.2.a 
where the pion vertex also includes the infinite set of counter terms needed 
to renormalize the bare result. The renormalized amplitude can be written 
as 
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Tabcd = 5 ab 5 cd A (s) + 5 ac 5 bd A (t) + 5 ad 5 bc A (u) 



(27) 



where 



with 



G r ( 8 ;vl) 



NF 2 i _ sGR ( s ^ T(s- u) 



mlG R {s-^) 
l-g^T(^) 



(4tt) 2 F 2 



log 



T(s;ii) = 



'i-Mi og 



1 - 



4mJ 



1 m 2 
- - log 



4 "4 i j 

s 



(29) 



The G R (s; ji) and J fi (s; /i) are the generating functions of the counterterm 
couplings and they have a well defined dependence on the \i scale, in order 
to ensure the \i independence of the whole amplitude. The isospin channels 
for the model here considered are defined as [14] : 



T = NA (s) + A (t) + A («) 
7\ = A(t) — A (u) 
T 2 = A(t)+A(u) 



and the partial waves are given by: 

1 n 



a u (s) 



64tt 



d (cos 6) Tj (s, cos 6) Pj (cos 6) 



(30) 



(31) 



As it is well known, the requirement of unitarity strongly constraints the 
possible behaviour of these partial waves. In particular, they should have a 
cut along the positive real axis from the threshold to infinity. The physical 
amplitudes are the values just on the cut, and, in addition, the condition of 
elastic unitarity: 



I ma! j = a\au\ 



(32) 



(where a = \Jl — 4m 2 / s) has to be satisfied in the physical region s = E 2 +ie 
and E 2 > 4m 2 , where E is the total center of mass energy. This equation is 
exact for energies below the next four pion threshold but even beyond that 
point it is approximately valid. 

In the level of approximation considered in this work (large iV limit and 
first order in m 2 /F 2 ), the dominant channel is the I — J — 0. As it can be 



8 



expected the elastic unitarity condition is not exact but it can be shown that 
one has 



From the partial waves it is possible to define in the standard way the phase 
shifts 5 as 



Now one can try to compare the results of our approximation with the exper- 
iment. Then one finds that the I = J = pion scattering phase shifts can be 
fitted with the simple choice G R (s; /i) — 1 and J R (s; //) = 1 for /j, — 775MeV, 
as it is shown in Fig. 3. Concerning the I = 0, J = 2 channel (which is the 
other relevant for the 77 — » 7r°7r° reaction), the prediction of the large N 
approximation is that it is suppressed at leading order in the 1/N expansion. 
In any case, one can use the above equations with the \i value fitted for the 
a o channel. The result, which is in fact a prediction, is also shown in Fig. 3 
and, as it can be seen, the agreement with the experimental data is pretty 
good. 

At the same level of approximation, the diagrams contributing to the 
77 — > 7r°7r° reaction can be found in Fig.2.b since it is not difficult to see 
that any other diagram is suppressed by extra 1/N or m\j F 2 powers. Thus 
the amplitude corresponding to this reaction can be written as 



AV = e 2 (2g,J(s) - I(hk 2 ),„ - I'(hk 2 )^) NA(s) + 0(1/N) + 0((-f ) 2 ) 




(33) 




(34) 




(35) 



where 




(36) 



The last two integrals have the same value, so that we can write 




fO(l/iV)+0((-^) 2 ) 
(37) 
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Note that this is the same integral that appears in the one-loop computation 
in eq.14. Therefore we also find, in the approximation here considered, the 
factorization property that we found there. In particular, eq.16 can be ap- 
plied again although now the lowest order 7r + 7r~ — > 7r°7r° cross-section has to 
be replaced by the cross-section obtained from eq.28, which is given by: 



^oA = I A ( g ) I' (38) 



The results obtained by using this formula, with the previously fitted value 
for pi in the I — J — elastic pion scattering channel, are given in Fig.4. 
For completeness we also show the results obtained when the J = 0, / = 2 
channel contribution is included by using the exact formula 

a (vr+vr- - vrV) = - ^ sin 2 (5 00 - 6 20 ) (39) 

together with eq.16 and the phase shifts computed in the large N approx- 
imation described above which are shown in Fig. 3. As it can be seen, the 
effect of the 1 = 2 channel is not very large, since it is not so strongly in- 
teracting as the 1 = 0. In any case our computations are compatible with 
the experimental data in both cases (with and without the introduction of 
the 1 = 2 channel) if one takes into account the large error bars. Thus, the 
large N approximation makes it possible to fit all the data shown in Fig. 3 
and Fig.4 with just one parameter, namely the scale // = 775MeV. 



5 Discussion 

In some sense the aim of this work is to clarify to what extent the above 
mentioned factorization relation between the 77 — > 7r°7r° and the 7r + 7r~ — > 
7r°7r° amplitudes occurs beyond the one-loop approximation. Note that, as 
it was stressed by Pennington and Morgan in [12], the computation of the 
exact amplitude will require the knowledge of the 7r + 7r~ — > 7r°7r° amplitude 
not only on-shell but also off-shell in order to make loops. In addition, the 
partial waves of the first reaction behave at threshold as s J//2 (s — 4m 2 ) J / 2 
whereas the second goes as (s — 4m 2 ) 3 and therefore the factorization formula 
cannot be correct. Moreover, from elastic unitarity (at the lowest order in 
the electromagnetic interactions), we have the relation 

lmj=l = aFfttu (40) 
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involving the 77 — > tt 7t° and the 7T7T partial wave (we have used the standard 
definitions which can be found for instance in [5]). In this reference it was 
shown how it is possible to use dispersion relations to improve the unitary 
behaviour of one- loop %PT, both for elastic pion scattering and 77 — > 7r°7r°. 
As a result, an excellent agreement with the data can be found in both cases. 
These dispersion relations have to be applied to each channel independently. 
Therefore, as the total amplitude is a linear combination, with well defined 
coefficients, of the different channel partial waves, we arrive again to the 
conclusion that the factorization relation cannot be hold in general. 

Therefore, for many different reasons, we known that the factorization 
formula which appears as a result of the large N approximation (as well as the 
original one-loop) cannot be true for the exact amplitudes. Then the natural 
question arises: Is there anything wrong with the large N approximation? 
and in that case: Why does it reproduce so well the experimental data? 

Of course the answer to the first question is no. The large N approxima- 
tion is a well defined expansion in the 1/N parameter and it is perfectly right. 
In particular the computations done by using this approximation satisfy the 
current constrains of any Quantum Field Theory like unitary, analyticity, 
crossing and so on, modulo higher order corrections in the expansion param- 
eter, i.e. 1/N (see eq.33), exactly as it happens in any other perturbative 
expansion. The reason why it gives rise to the factorization formula is the 
following: according to eq.27 the A function, which in principle depends on 
s,t and u, is only a function of s at leading order in the 1/N expansion. In 
addition, from eq.30 we see that, strictly speaking, only the 1 = channel 
contributes in this approximation. In other words we can say that to leading 
order in the large N expansion the only non-zero channel is I — J — 0. Then 
the 7r + 7r~ — > 7r°7r° amplitude only depends on s and only has an isoscalar 
component. This explains why all the previous mentioned objections to the 
factorization formula do not apply to this case as it happens in the general 
case where more than one channel contributes to the reaction. Therefore 
there is no inconsistency in the factorization formula obtained in the large 
N limit. However, this formula should be corrected by higher order compu- 
tations (both in 1/N and in m 2 /F 2 ) which will also include other channels 
different from the I — J — 0. In that case the factorization formula is not 
justified any more. 

Of course one can ask how good the large N approximation to leading 
order in 1/N and in m 2 /F 2 can be, since it only includes the I — J — 
channel contribution. The answer can be found in Fig.4. There we see that 
the agreement with the experimental data is good, and in fact it considerably 
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improves the one-loop computation (the two-loop computation in [5] also 
reproduces well the data but depends on more parameters) 

6 Conclusions 

The main results of our work are the following: 

In order to properly define the large N expansion for two flavor %PT we 
have parametrized the coset space S* 3 — > S N by using standard coordinates 
instead of exponential (chiral) coordinates. We have reobtained the well 
known one-loop results for the 77 — > nn reactions with this parametrization 
in a much simpler way, thus illustrating, with a very non-trivial example, the 
S matrix invariance under reparametrizations of the NLSM. 

We have computed the 77 — > 7r°7r° scattering amplitude to leading order 
in the large N approximation up to order ml/F 2 . In this approximation 
this amplitude turns out to be proportional to that of it + it~ — > 7r°7r° (the 
factorization formula) according to Im's conjecture (but in principle only for 
small ml/F 2 ). 

Using a one-parameter fit of the elastic pion scattering in the I = J = 
channel obtained in the same approximation we reproduce, as a prediction, 
both the I = 2, J = pion phase shift and the 77 — > tt°-k° cross-section 
experimental data (see Fig. 1 and 2). 

From very general considerations it can be shown that the factorization 
formula cannot be exact and this was the basis for the Morgan and Penning- 
ton criticism to Im's results. However Im's result is correct in the framework 
of the large N approximation (for small ml/F 2 ) since only the I — J — 
channel contributes to the leading order in the 1/N expansion. As a matter 
of fact, as the 77 — > 7r°7r° is dominated by this channel, the factorization 
formula yields a good approximation as it was suggested by Donoghue et al 
in [3]. This would not be the case if, for example, the p resonance played an 
important role in the reaction, since it cannot be reproduced by the leading 
order large N approximation [11]. 

In summary, we consider that the large N approach to is a very 

interesting complementary alternative to the standard one-loop computations 
which can be useful in many cases in order to describe the relevant physics 
of the considered phenomena. This is particularly true in those cases where 
the I = J = is dominant as it happens in the Higgs physics. 
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Figure Captions 

Fig. 1: Diagrams contributing to a) the 77 — > 7r°7r° reaction at one- loop 
level, b) the tree level 77 — > 7r + 7r~, c) the one-loop 77 — > 7r + 7r~ and d) the 
one-loop pion wave-function renormalization. 

Fig. 2: a) Diagrams contributing to the 7r — 7r elastic scattering at leading 
order in the 1/N expansion up to order m 2 /F 2 . The large black dot repre- 
sents the pion coupling proportional to m 2 . b) Diagrams contributing to the 
77 — > 7T 7r° reaction to the same level of approximation. The black square 
represents the addition of all the diagrams appearing in a). 

Fig 3: Phase shifts for J = elastic pion scattering. The dashed line 
represents the fit done with the leading order of the 1/N expansion up to the 
order m 2 /F 2 for the 1 = channel and the prediction for the 1 = 2 channel. 
The continuous lines are the phase shifts as obtained from the standard one- 
loop X-PT (see [11] for the details). The experimental data corresponds to: 
A ref.[15], O ref.[16],D ref.[17],0 ref.[18], y ref.[19],* ref.[20],x ref.[21],« 
ref.[22]. 

Fig 4: Cross section for the process 77 — > 7r°7r° with | cosO \< 0.8. 
The continuous line represents the prediction of leading order of the 1/N 
expansion up to the order m 2 /F 2 . The dashed line represents the same but 
including the contribution of the 1 = 2 channel. The experimental data 
comes from Cristal Ball [1]. 
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